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Abstract 

We investigate the evolution of a given eigenvector of a symmetric (de- 
terministic or random) matrix under the addition of a matrix in the Gaus- 
sian orthogonal ensemble. We quantify the overlap between this single vector 
with the eigenvectors of the initial matrix and identify precisely a "Cauchy- 
flight" regime. In particular, we compute the local density of this vector in 
the eigenvalues space of the initial matrix. Our results are obtained in a non 
perturbative setting and are derived using the idea of Ledoit and Peche in 
[llj. Finally, we revisit our former results on the eigenspace dynamics giving 
a robust derivation of a result obtained in pQ in a semi-perturbative regime. 

1 Introduction 

The dynamics of eigenvalues induced by the addition of free random matrices in the 
Gaussian orthogonal ensemble has been first studied by Dyson in his 1962 paper [8] . 
The movement of the eigenvalues is characterized in terms of a stochastic differential 
system, the so called Dyson Brownian motion. The eigenvalues evolve as particles 
of a Coulomb gas with electrostatic repulsion, confined in a quadratic potential and 
subject to a thermal noise. In the limit of large dimensions, the evolution of the 
spectral density has also been studied in [T3] (see also [51 E] and [2] for related 
models). 

For the eigenvectors, their evolution in finite dimension is also given by a stochas- 
tic differential system which depend on the non colliding trajectories of the eigen- 
values (see |3]). In this paper, we are interested in quantifying the evolution of the 
eigenvectors in the limit of large dimension. Our approach uses the idea of [UJ who 



introduced a very interesting quantity (see Eq. (5.1) below) for the study of eigen- 
vectors. This enables us to compute the local density of a given state (eigenvector) 
of the matrix after the addition of the free Gaussian matrix, in the eigenvalues space 
of the initial matrix. 

The paper is organized as follows. In section |2j we define the model and give the 
main notations. We enunciate our main result (Theorem 3.2) in section [3] concerning 



the convergence of the quantity (5.1) introduced in [11J. In the following section H 



we are interested in the local density of states under free addition. We find that in a 
particular regime, the eigenvalue dynamics can be precisely described as a "Cauchy 
flight" . We also check numerically our results in the case of a initial random matrix 
in the Gaussian orthogonal ensemble. We then revisit in section 4A the main result 
of pQ on the dynamics of eigenspace under free addition, and prove that it is indeed 
exact beyond the perturbative regime. Finally, section [5] is devoted to the proof of 
Theorem 13.21 
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like to thank also Sandrine Peche for interesting discussions. 



2 Definition and main notations 

Let A be a symmetric deterministic N x N matrix. The matrix A is diagonalizable 
in an orthonormal basis of M. N and its eigenvalues will be denoted in increasing 
order as a± ^ a 2 ^ • • • ^ ajy- We will suppose that the dj are allocated smoothly 
as = a(-h) where a(-) is a "nice" differentiable and strictly increasing function 
defined on the interval [0; 1]. 

Let (H t ) t ^ o be a symmetric Brownian motion, i.e. a symmetric diffusive matrix 
process constructed from a family of independent real Brownian motions Bij(t), 
l^i^j^N&s follows 

H t(ii) M Bi ' (t) m<i ' 

(y$fls(t) »<->. 

Note that the process H t is defined such that it is rotationally invariant at all time, 
in the sense that for all O in the orthogonal group On, the conjugate matrix OH t O^ 
has the same law as the matrix H t . 

Now we define the matrix M t such that for all t ^ 

M t = A + H t . (2.1) 

The eigenvalues of the matrix M t will be denoted in increasing order as ^ • • • ^ \\ 
The aim of this note is to quantify the relationship between the eigenvectors of 
the matrix M t with the eigenvectors of the initial matrix M = A. In particular, 
we consider one given eigenvector of the matrix M t denoted as \ipj) and we want 
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to compute, in the large limit of large dimension N, the overlaps of this vector 
with the eigenvectors of A denoted in the following as \(f>i), ■ ■ ■ , \4>n) (which form 
an orthonormal basis). As the matrix H t is rotationally invariant for all t, we can 
suppose with no loss of generality that the matrix A is diagonal. 



3 Results 

If t is very small, perturbation theory for the eigenvectors of the matrix M t defined 



in (2.1) gives an approximation of the vector \ip\) in terms of the small matrix H t 
and of the eigenvectors of the matrix A. More precisely, if t is small compared to 
the level spacings of the matrix A, perturbation theory leads to (see appendix[A|for 
a reminder of this computation) 

/ i r (^wV, u viMu\^i f \ 

^ = I 1 " 2 \ (a, - a,) 2 j + \ ~^a~ ^ + N ' (5U) 
where \e t ) is a vector such that {(pi\e t ) = o(t) and (<pj\e t ) = o(\/i) for j ^ i. This 



expansion (3.1 ) can be rewritten, taking a further expectation, in terms of the overlap 



between and \<pj) for i ^ j as 
and for i = j, 



The computations leading to equation (3.1) are a priori only valid for values of 



t much smaller than the spacings between the consecutive eigenvalues a, L — a^i, 



so that the correction terms in the perturbation equation (3.1) remains very small 
compared to the leading term of order 1. Note that in our setting and for large N, 
a.i — cij_i is of order 1/N, so that t has to be much smaller than 1/N. 

The purpose at stake here is to handle the case of larger values of t in the large 
limit. For this eigenvector dynamics problem, there are in fact three regime for 
the values of t and A^: 

• The first regime is the perturbative regime where t 1/A^, for which the 
perturbation theory applies to give the approximation ( 3.1[ ). This regime is 



well known and has been studied in great details in random matrix theory and 
in the context of quantum mechanics. 

The second regime is what we will call semi-perturbative and concerns small 
values of t compared to 1 but not necessarily small compared to the level 
spacings of order 1/N . It includes values such that 1/N <Ct < 1. 
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The third regime is non perturbative: N is very large while t is of order 1. 



The question we ask is: How do formulas (3.2) and (3.3) get modified in the second 



and third regime in the large N limit ? Because the family {cj)j, 1 ^ j ^ N} forms 
an orthonormal basis of Mr, we have 



N 



E^%> 2 = L 



=i 



We will therefore investigate the convergence in the large N limit of -/VE[(?/^|0j) 2 ] 
for i j. Those quantities are in fact related to the mean local density of state \ip\) 
which is the probability measure Ui defined on the ^-eigenvalues space as 

N 
3=1 

where 8 (x) is the Dirac measure in x. In other words, the aim of the note is to 
compute the local density vi of state \ipf) in the large N limit. 

The interesting quantity for our purpose is the bivariate cumulative distribution 
function $ associated to the weights NK[(xjjj \4>j) 2 } defined for A, a G M by 



N 

hi 



1 N 



{A* < A} ^{a, < a} 



Note that this function $ has indeed the properties of a bivariate cumulative dis- 
tribution function since it is right continuous with left-hand limits and is nonde- 
creasing in each of its argument and that it satisfies lim $(A,a) = and 

A— >— oo,a— >— oo 

lim $(A, a) = 1. 

A— >+oo,a— >+oo 

Before presenting our results on the convergence in the large iV limit of the 
bivariate cumulative distribution $jv(A, a) which will directly lead us to asymptotic 
estimates for the overlaps iVE[(-?/^| (/>,,• ) 2 ] for i ^ j and for the local density of states, 
let us first recall a result on the convergence of the empirical eigenvalue distribution 



of the matrix M t defined in (2.1) due to Shlyakhtenko [H] (see also [X2| [19] for 
similar results). Recall the definition of the Stieltjes transform G^(z) of a probability 
measure fi on R defined for all z G C \ K as G^(z) = J R . 

Theorem 3.1 (Shlyakhtenko, [H]). The empirical eigenvalue distribution of the 
matrix M t converges weakly almost surely when N — > oo to a probability measure \i t 
whose Stieltjes transform satisfies 

G m( z ) = / G t (z,x)dx 



o 
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where G t {z,x) is defined uniquely for any z G C \ R by the fact that it is analytic 
on C \ R ; maps i/ie upper half plane C + mto £/ie lower one C _ , and satisfies the 
relation for all z G C \ R ana 1 x G [0; 1] 

G t (z,x) = — • (3.4) 

a(x) - z -tj Q G t (z,y)dy 

Furthermore, setting Rf(z) = (M t — zl)" 1 , the complex random measure on the 
interval [0; 1] defined as 

1 N 

8=1 

converges weakly almost surely to the complex measure with density G(z,x)dx. 

The proof of this theorem goes through the Schur complement formula, the fixed 
point theorem and classical properties of analytical functions. 

Finally, using the work of Biane [Bj, we know that the limiting spectral distri- 
bution p t (dX) admits a smooth density pt(A) with respect to Lebesgue measure. We 
further stress the convergence of the Stieltjes transform near the real axis: 

lim Gat(X + it]) = H Pt (X) + inp t (X) 

??->0+ 

where H pt is the Hilbert transform of the probability density p t . 

We are now ready to state our main result on the convergence of the bivariate 
cumulative distribution $jv(A,a). 

Theorem 3.2. When N — > +oo, we have the following convergence 

f x f a dx t 

$jv(A,a) — ><5>{X,a)= / d£p t (£) 



a'(a-i(z)) (x - C ~tH pt (0) 2 + tVp t (0 



We will see in the following section that Theorem |3.2 permits to compute the 



asymptotic overlaps A^E[(^*|0j) 2 ] of any eigenvector of the matrix M t at time t 
with the eigenvectors of the initial matrix Mq = A for any time t. A former heuristic 
attempt to compute those overlaps already appeared in [17], with a different result 



(see below). Theorem 3.2 also permits to compute the mean local density U{ in the 
^-eigenvalues space of the state \ipf) at time t. In addition we will explain how 



Theorem 3.2 enables to extend the domain of validity of our former result on the 



eigenspace dynamics under free addition obtained in [TJ. 



Theorem 3.2 is similar to [TT] Theorem 1.3]. In [TT], the authors investigate the 
relationship between the eigenvectors of the population covariance matrix with those 
of the empirical (or sample) covariance matrix. 



One remarkable feature of Theorem 3.2 is that it quantifies the relationship 



between the eigenvectors of the initial matrix Mq = A and the eigenvectors of the 
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matrix M t , even in the non perturbative third regime described above (i.e. it holds 
for any time t > 0). In [TT], the sample covariance matrix is an additive perturbation 
of the population covariance matrix. The authors of [TT] thus work in the second 
semi-perturbative regime described above (note that there is not a third regime in 
their case though). 

Our proof of Theorem of 3.2 uses the idea used in [11] to quantify the relationship 
between sample and population eigenvectors. 



4 Applications of Theorem 3.2 



The present section contains physical arguments and is not fully rigorous mathe- 
matically. 



4.1 Perturbation theory revisited 



Theorem 3.2 enables to compute the overlaps of the vector \ipf) with the initial 
(t = 0) eigenvectors \<f>j) for j ^ i of the matrix Mq = A. We can deduce from 
Theorem |3.2 that, for N — > oo, we have 



(4.1) 



As mentioned before, formula (4.1) is valid for all values of t and N in the third 



regime described above (i.e. for large N and any t > 0). 

For f<l and for eigenvalues a i: aj separated by a macroscopic spacing (i.e. such 
that a 

as 



<Zj does not vanish for large N), we have A* ~ a, and formula (4.1 ) rewrites 

t 1 



+ <]v>- 



(4.2) 



For such i,j, (4.2) extends the perturbation equation 3.2| (which is valid only for 



i < 1/ N) to values of t much smaller than 1 but not necessarily negligible compared 
to 1/N. 



4.2 Local density of state 



In other words, Theorem |3.2 gives the local density of states in the large N limit. 
The limiting mean local density of states |^*) in a is 

Vl{a) = a>(a-i(a)) (a - X t - tH pt {\)) 2 + t^p t {\Y ' (43) 



This last formula (4.3) defines a probability density on K although it is not trivial to 



check that its integral over a is indeed 1 (the probability measure p t and its Hilbert 
transform depend on the function a(x) so that ^ is indeed a probability density 
function). 
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4.3 The case of Dyson Brownian motion 

In this section, we apply our results to the particular case in which the initial matrix 
A is a random matrix of the Gaussian orthogonal ensemble (GOE). As before, we 
will work in the basis of the eigenvectors of A, so that the matrix A is diagonal 
with eigenvalues in increasing order: A = Diag(ai, ■ ■ • , a^) with a% < ■ • ■ < a^. The 
spectrum of A is random in this case, and the assumption we made to derive our 
results on the smooth allocation of its spectrum is not satisfied here. We will check 
numerically that our results hold also in this case though. 

Using Wigner's semi circle law, we can conjecture that the spectrum of A should 
converge to the continuous and strictly function a : [0; 1] — > [—2; 2] such that, for all 
XG[0;1], 

I r a ( x ) 

— / a/4 - t 2 dt = x. 
2n J_ 2 

The matrix M t is also a GOE matrix and its limiting empirical eigenvalue density 
is given by the semi-circle law 

"' (A) = MTT7)^ 4(1 + t) - A2 - (44) 

It is also known (see e.g. [3]) that the associated Stieltjes transform G^iz) satisfies 

lim Gfj, t (X + irj) = - X + mp t {X) ■ 
ii->o+ 2(1 + t) 



Hence, Theorem |3.2| gives the following expression for the asymptotic overlaps for 
«V 3- 

NEMMj)*] = 7 i * 777—1^ + o(l) ■ (4.5) 

We have checked those asymptotic expressions, using numerical simulations, in the 
present context of GOE matrices of dimension iV = 400 with t — 1. The agreement 
is excellent: see Fig. [T]and[2j 

In other words, the local density Vj of state \if)l) in the limit of large matrices is 
the probability measure with density in the A-eigenvalues space 



. . \/4 - a 2 t 



When t — > +oo, the information about the initial state is finally lost 

lim lim A^%> 2 ] = 1. 

t— >-+oo N— >+oo 

The vector \ipf) has uniform overlaps with the initial eigenvectors \<j)j) as expected. 
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In the limit tcl, Eq. (4.5) can be simplified to read: 



(oj - A 



t\2 



+ t 2 



(4.6) 



which describes a "Cauchy flight" in eigenvalue space. This makes more precise a 
statement made in [TTt H~8| [I] in the context of an extreme nonadiabatic evolution of 
a quantum system: the energy is not diffusive but rather performs a Cauchy Flight. 

In fact, if the evolution of the system is such that the elements of the random 
GOE matrix M t have a fixed variance, the z'-th eigenvalue of M t is expected to be time 
independent in the large N limit, i.e. A* ~ a,. In this case, Eq. (4.6) corresponds (up 



to simple modifications) to Eq. (4.11) of |17j . with the correspondence AE = a, — a*. 
However, the correspondence for longer "times" t must take into account that with 
our normalization, the semi-circle spectrum itself broadens with time, as given by 



Eq. (4.4). 



o 



in 
o 




Figure 1: The black curve is a plot of N~E[(ipj\(f)j) 2 ], computed empirically with 1000 
samples, as a function of the eigenvalues dj corresponding to \4>j), for N = 400, t — 1 
and % = 200. The 200 th eigenvalue of A is approximately equal to 0, it is therefore 
natural to observe the highest value of this curve at this point. The red curve is the 



theoretical prediction displayed in Eq. (4.5). 
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Figure 2: The black curve is a plot of iV E [(ipj \<j)j) 2 }, computed empirically with 
1000 samples, as a function of the eigenvalues Oj corresponding to \<fij), for N = 400, 
t — 1 and i = 320. The 320 th eigenvalue of A is approximately equal to 0.983, it is 
therefore natural to observe the highest value of this curve near this point. The red 
curve is the theoretical prediction displayed in Eq. (4.5). 
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4.4 Eigenspace stability 

In PQ, we investigated the stability of eigenspaces associated to a GOE matrix A 
when a small GOE matrix H t is added. Let us briefly recall the context and main 
notations of [I]. 

Our idea was to study, in the large N limit, the stability of a whole subspace 
Vq (instead of a single eigenvector as above) spanned by a set of consecutive initial 
eigenvectors \4>k) associated to eigenvalues contained in a certain interval py_; 7+] 
of the Wigner semicircle support [—2; 2]. We then asked the following question: How 
should one choose a "larger" subspace V( spanned by a subset of eigenvectors j^) at 
time t which would contain the initial subspace Vq up to a small error ? To answer 
this question, we introduced a margin of width 6 and the subspace V{ generated 
by the set of eigenvectors \ipl) associated to eigenvalues A| lying in the interval 
[7_ — 5; 7+ + 5] . We then considered the rectangular matrix of overlaps Gt with 
entries 

G t (ij) :=(^%>- 
In this setting, the matrix Gt has dimensions Q x P with 

P = N po, Q = N p t , 

J 7_ J , y——5 

where po is the Wigner semicircle eigenvalues density of the initial matrix A. The 
labels i and j and the vectors \if)f) and \4>j) are respectively indexed by the eigenvalues 
(in increasing order) A- and aj. 

The P nonzero singular values 1 ^ Si ^ S2 ^ ••• ^sp^Oof the matrix 
Gt contain full information about the overlap between the two spaces Vq and V{. 
For example, the largest singular value S\ indicates that there is a certain linear 
combination of the Q eigenvectors at time t that has a scalar product s\ with a 
certain linear combination of theP initial eigenvectors. If Sp — 1, then the initial 
subspace is entirely spanned by the perturbed subspace. If on the contrary Si <C 1, 
it means that the initial and perturbed eigenspaces are nearly orthogonal to one 
another since even the largest possible overlap between any linear combination of 
the original and perturbed eigenvectors is very small. A good measure of what can 
be called an overlap distance D(Vq, V r /)between the two spaces Vq and V\ is provided 
by the average of the logarithm of the singular values: 

1 p 

D(Vo,Vf) = --^ln( Sfc ) . 

k=l 

This overlap distance D and the overlap matrix Gt already appeared in the literature 
on the "Anderson orthogonality catastrophe" (see e.g. jHHUj). 
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Using perturbation theory, we showed in [T] that this overlap distance D(Vq, V*), 
ioT N ^> 1 and in Nt <C 1, converges to: 

E[D(V , VI)] = -±- [ J+ dx [ dy P f )p ^ + o(t) . (4.7) 

We conjectured that this formula also holds in the semi-perturbative regime N 3> 1 
and t C 1 but Nt ^> 1. This conjecture was supported by convincing numerical 
evidence pQ . Unfortunately we were unable at the time to find analytical arguments 
to sustain our claim in this semi-perturbative regime. 

We will use our new results obtained here to fill this gap and prove that Eq. 
(4.7) is indeed correct in the whole regime for i < 1 (independent of N) in the limit 
N -> +oo. 

Let us first remark that D(V 0l Vf) = - ln(det(GjG t ))/(2P) where G\ is the 
Hermitian conjugate of Gt- We thus compute the entries of the matrix G t Gt- For all 
di G [7_;7_|_], we have 

(G\G t h = WM' 
A' fc e[ 7 --5;7++<5] 

= <$l*> 2 + E M**)' 

= i-E<^> 2 + E <^> 2 

k^i fc^j:A t fc e[7--'5;7++5] 

= i- E wm 2 - 



Using the previous result (4.5) obtained in the previous subsection, we see that for 
values of t <C 1 and for large N, we have 

a fe 0[7_-5; 7+ +5] 

-at^oo 1 - t I dy Poiy \ 2 + o(t) (4.8) 

Jy^[7_-<5, 7+ +<5] l a i — V) 

where we have used the fact that for small values of t, the eigenvalues satisfy ~ 
A° = afc for the last line. 

The non diagonal elements, i.e. indexed by a, ^ £ [7_;7 + ], can also be 
computed as 

{G\G t ) l3 = WMWM 

A' fc e[ 7 --<5;7++<5] 

= - E WMMMj) 

A^[ 7 _-<5; 7+ +«5] 
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where, in the second line, we have used the orthogonality of and \<pj) which 
implies that Y^ki^kl&i) (^kl&j) = ®- The expectations of those terms can thus be 
estimated for small values of t and large N, via the Cauchy-Schwarz inequality, as 



E[(G\G 



t 



/2 



N 



E 



1 



-5;7++<5] 

dy 



(at - a k )(aj - a k ] 

po(y) 



+ o{-)+o{t) 



y ^_-5, J+ +6] K - y)( a j - v) 



+ 0(t) . 



(4.9) 
(4.10) 



Note that in both equations ( |4.8 ) and (4.10), we have a^aj E [7-; 7+] and a k £ 
[7_ — S;j + + 5), so that a; (resp. a,) and remain at macroscopic distance > 5 
and formula ( 4.2[ ) applies. Besides the integrals in (4.8) and (4.10) are perfectly well 
defined due to the introduction of the margin 5 > 0. 

Thus, for t <C 1 and large N, the determinant of G\Gt can be approximated, in 
the large N limit, to leading order in t as the product of the diagonal terms (the 
other contribution are negligible compared to t for small t). We thus have, doing a 
further linearization in the limit of small t, 



-— Epn(det(G{Gt))] 



t 



E , 

a i e[7-;7+r I/ * [7 — 5 ' 7++5] 

t /" 7+ 



A i / 1 \ , f+\ 



n+ 
J 7 _ 



Po 



dx 



dy 



j/£[7_-<5,7 + +(5] 



Po{x)po{y) 

(x - y) 2 



o(t). 



This is our proof that (4.7) is valid in the second semi-perturbative regime. 



The reader may wonder how to extend formula (4.7) in the non perturbative 
regime, i.e. for arbitrary values of t. This question is clearly more difficult as one 
would need to understand the convergence of the non diagonal terms of the matrix 
Gt f Gt in the large N limit which are no longer negligible in the determinant 
expansion. 



5 Proof of Theorem 13.2 

Following the idea of [H] , we introduce the following quantity, defined for z G C \ R, 

as 

N N 
i=l 1 j=l 

= ^Tr ((M t - zI)- x g{A)) (5.1) 
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where g is a real valued bounded function on R. By convention, g(A) is the diagonal 
matrix Diag(</(ai), 5-(a 2 ), • • • ,g(a N )). 

The interesting feature of Q g N {z) is that, by the Stieltjes inversion formula, we 
have 

$ N (\,a) = lim - / Im[e 9 N (£ + ir])}d£ 

»7->0 + 7T 

for the particular choice g(x) = 1{ X ^ a }- 

Thus, the problem is reduced to the study of the convergence of Q 9 N {z) when 
iV — > oo. It is plain to deduce Theorem 3.2 from the following lemma. 

Lemma 5.1. Let g be a real valued bounded function on R. Then, as N — )■ +oo, we 
have the following convergence 



e 9 M 



o(a(a;)) 



a(a;) — z — tG^iz 



-dx 



where G^ t (z) is the Stieltjes transform of the limiting spectral distribution fi t of the 
matrix M t . 



Proof of Lemma 5. 1 



Using equation (5.1) and the definition of the matrix (z) = (M t — zl) , it is 
straightforward to check that 



N 



i=l 



I 

N' 



Now, using Theorem 3.1 of Shlyakhtenko (see [H], and also [T2"| [TU]). we know that 
the complex- valued measure af defined in (3.5) converges weakly to G t {z,x) dx. 
Therefore, as N — > oo, 



g(a(x))G t (z, x)dx 



g{a{x)) 



a(x) — z — tG m (z 



-dx . 



using the fixed point equation (3.4) satisfied by Gt(z,x). The lemma is proved. □ 



A A reminder of perturbation theory 



We can rewrite (2.1) as 

M t = A + Vt Hi 

where Hi is a GOE matrix whose upper diagonal entries are independent Gaussian 
variables of variance 2/N on the diagonal and variance 1/N off the diagonal. 
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To second (respectively first) order in yt, we expect the following perturbative 
equations for the eigenvalue A* and the associated eigenvector of the matrix M t 

\ t i = a i + y/ia i + t/3 i + o(t) , (A.l) 
|$) = (1 - 7i t) |&> + Vi J2 7j \<P 3 ) + \e t ) . (A.2) 

where \e t ) is a vector such that {4>i\s t ) = o(t) and (0j|£t) = o{\fi) for j 7^ i. 

The aim here is to compute the terms ati, fli and 7, for j i. We start writing 
the eigenvalue equation for the matrix M t 

M t |#)=A*|$) 



which can be rewritten plugging equations (A.l) and (A.2) as 

y^7j|0j) + t# +ai|e t ) + |ei) , 

where \e t ) is a vector such that {4>i\st) = o(t 3 ^ 2 ) and (4>j\e t ) = o(t) for j ^ i. 

Projecting this equation on |0jrjand identifying the leading term in y/i on both 
sides, we find 

on = (^i\Hx\(f)i) ■ 

Then, projecting on \<pj) j ^ i and identifying the leading term in \/t on both sides, 
we find 

\<l>i\ H Mi) 



di — CLj 



Projecting again on \<pi) but identifying the term in t, we obtain 

The coefficient ji is computed with the additional normalization constraint = 
1 to first order in t. 

This gives the perturbation equations for the eigenvalue A- and the associated 
eigenvector Note that this perturbative computation is valid only if the terms 
t^ijtji and y/i'Jj remain very small compared to 1, corresponding to small values 
of t compared to the level spacings a,i — aj (which are of order 1/N in our setting 
above). 

1 Note that the vectors \cj>k) for k = 1, • • • ,N form an orthonormal family. 
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